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Abstract 

This  report  deals  with  problems  of  transport  properties  in  gas  dynamic  flows  at  high  pressures  and  densi- 
ties and,  in  particular,  with  transport  properties  of  Hydrogen,  Nitrogen,  Water  and  Oxygen.  The  behavior 
of  these  substances  at  high  pressures  and/or  low  temperatures  is  quite  different  from  perfect  gas  behavior. 
Especially  near  its  critical  points  or  lines  of  phase  changes,  thermodynamic  properties,  equations  of  state 
and  transport  properties  as  well  exhibit  proper  characters  while  the  specified  conditions  are  taking  place. 
An  overview  of  models  in  use  to  predict  the  transport  properties  for  the  species  given  above  in  various 
phase  states  is  carried  out  and  different  theoretical  approaches  are  investigated.  A generalization  of  the 
models  requires  the  knowledge  of  all  these  properties  in  the  widest  possible  range  of  temperatures  and 
pressures.  Therefore,  the  inclusion  of  limiting  states'  data  (dilute  gas,  saturation  and  melting  regime)  into  a 
multi-property  analysis  is  essential,  especially  while  extending  these  correlations  beyond  the  range  of  ex- 
perimental data. 

A novel  method  for  the  prediction  of  viscosity  and  thermal  conductivity  on  the  basis  of  an  unique  equation 
for  gas  and  liquid  is  presented  which  bases  on  the  model  of  effective  kinetic  diameter  of  hard  spheres. 
Furthermore,  the  kinetic  effects  on  the  droplet  surface  under  high  ambient  temperature  are  considered. 
Temperature  and  concentration  jumps  were  used  to  describe  non-equilibrium  boundary  conditions  on  an 
evaporating  liquid  oxygen  (LOX)  droplet.  The  influence  of  calculation  errors  on  the  values  of  these  jumps 
is  evaluated. 

Introduction 

The  properties  of  matters  in  various  phase  states  may  be  divided  into  two  kinds,  equilibrium  properties 
and  transport  properties.  A great  progress  has  been  made  in  the  recent  years  in  studying  equilibrium  prop- 
erties, although  inconsistency  between  calculated  properties  is  often  a problem  in  industrial  thermo- 
physical properties  simulations  because  typically  different  prediction  methods  are  used  for  different  phase 
states.  Hence,  it  may  happen  that  equations  for  liquids  are  used  for  the  prediction  of  properties  in  a region 
where  the  vapor  pressure  equation  indicates  a gas  phase  should  be  present.  Therefore,  a consistent  method 
for  the  prediction  of  all  the  properties  based  on  a single  set  of  fundamental  parameters  is  the  most  prefer- 
able. A generalized  equation  of  state  can  be  used  accurately  to  predict  equilibrium  properties  of  fluids  and 
dense  gases. 

Furthermore,  transport  properties  of  matters  in  different  thermo-physical  states  are  also  important  features 
required  in  various  engineering  design  problems  such  as  simulations  of  viscous  flows  through  channels 
and  combustion  chambers  of  various  technical  devices  such  as  flows  in  rocket  engines,  chemical  reactors 
or  shock  tubes.  The  governing  equations  for  these  gasdynamic  systems  are  the  Navier-Stokes  equations 
with  mass  diffusion,  heat  flux,  pressure  tensor  expressions.  All  the  coefficients  in  these  equations  have  to 
be  known  precisely  since  they  have  a very  strong  influence  on  the  accuracy  and  the  consistency  of  the 
simulations. 

The  object  of  our  investigations  presented  here  is  (1),  to  provide  adequate  thermo-physical  property  data, 
mainly  transport  properties,  covering  the  largest  possible  range  of  temperatures  and  densities,  and  (2),  to 
demonstrate  a model  of  vaporization  accounting  for  thermal  non-equilibrium  boundary  conditions  on  the 
surface  of  a liquid  droplet. 
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Thermodynamic  properties 

Under  various  boundary  conditions  different  physical  phases  such  as  dense  gases  and  liquids  usually  co- 
exist, i.e.  the  flow  of  a mixture  of  cold  gaseous  hydrogen,  superheated  steam  and  liquid  oxygen  in  a cryo- 
genic liquid  rocket  engine.  Although  mixtures  of  layers  of  liquids  and  dense  gases  may  be  treated  as  a 
continuum,  the  properties  of  substances  under  high  pressure  and  low  temperature  are  quite  different  from 
those  of  a perfect  gas.  Hence,  thermodynamic  properties,  equations  of  state  and  transport  properties  as 
well  have  their  particular  characteristics  at  specific  conditions. 

The  conventional  interpretation  of  the  specific  characteristics  of  thermo-physical  properties  of  matters  un- 
der different  boundary  conditions  of  pressures  and  temperatures  generally  makes  use  of  well-known  phase 
diagram  of  states.  Figures  1 and  2 shown  typical  phase  diagrams  for  N2  and  H2. 


mc 

Fig.  1 : Phase  diagram  of  N2 
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Fig.  2:  Phase  diagram  of  H2 


All  three  different  physical  phases:  gas,  liquid  and  solid  states  are  realized  under  different  pressures  as  the 
temperature  varies.  The  triple  point  (TP)  is  the  point  of  co-existence  of  all  three  phases.  The  melting  line 
(ML)  is  the  boundary  line  between  the  liquid  and  the  solid  state.  The  sublimation  line  separates  the  co- 
existence of  gas  and  solid  states.  The  saturation  line  (SL)  is  the  boundary  line  between  gas  and  liquid, 
along  which  specific  volumes  of  gas  and  liquid  co-exist  at  given  temperatures  and  pressures.  The  critical 
point  (CP)  is  characterized  by  the  critical  values  of  temperature  T,  pressure  P and  density  pc.  The  CP  is 
the  final  point  on  the  SL,  where  the  density  (or  the  specific  volume)  of  a dense  gas  becomes  equal  to  the 
density  (or  the  specific  volume)  of  a liquid.  Above  this  point,  phase  transitions  don’t  happen  anymore.  The 
distinctions  between  gas  and  liquid  have  all  disappeared  above  CP,  and  the  specific  dense  gas  behavior  is 
close  to  that  of  a liquid.  In  other  words,  it  is  impossible  to  distinguish  between  gas  or  liquid  exists.  The 
choice  of  boundary  between  gas  and  liquid  has  become  a matter  of  convention  and  usually  the  critical  iso- 
chore  is  used. 
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132.62 
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273. 
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221.39 
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Table  1 : Critical  parameters  of  substances  of  interest 
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A brief  summary  of  the  critical  parameters  for  some  of  the  substances  of  interest  is  given  in  Table  1,  much 
more  details  can  be  found  in  [1].  Although  the  critical  parameters  for  various  matters  differ  widely  the 
particular  behavior  of  thermodynamic  (transport)  properties  near  critical  areas  is  quite  similar. 

Various  papers  have  been  devoted  to  the  investigation  of  equations  of  state  for  simple  fluid-gas  substances 
such  as  H2,  O2,  N2,  or  H20  in  wide  ranges  of  temperature  and  pressure.  The  outcome  of  both  experimental 
and  theoretical  treatments  made  it  possible  to  work  out  accurate  correlations  for  basic  tables  of  thermody- 
namic properties,  including  specific  volume  (V),  enthalpy  (/?),  entropy  (5),  and  specific  thermal  capacity 
( Cp ) in  terms  of  temperature  T and  pressure  P,  as  unique  reference  sources  of  data  [1][2][3].  As  an  exam- 
ple, the  specific  behavior  of  factor  compressibility  Z for  Oxygen  is  shown  in  Figure  3,  with  p the  pressure 
and  V the  specific  volume. 


pV 

Fig.  3:  Compressibility  factor  Z = - — 
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Fig.  4:  Thermal  Capacity  near  critical  point 


Near  CP,  a steep  increase  of  the  compressibility  factor  Z and  an  anomaly  of  the  thermal  capacity  can  be 
observed.  Fluctuations  in  density,  molar  fractions  and  other  physical  values  cause  a delay  of  diffusion,  a 
growth  of  viscosity  and  thermal  conductivity  and  give  rise  to  other  problems.  The  thermodynamic  proper- 
ties near  CP  were  investigated  for  H2,  02,  N2,  and  H20  thoroughly  and  may  be  found  in  [1]  or  [3].  An  ex- 
ample for  the  anomalous  behavior  of  the  specific  thermal  capacity  of  02  near  critical  region  is  shown  in 
Figure  4. 

Transport  Properties 

For  dense  gases,  viscosity,  thermal  conductivity  and  thermal  diffusion  as  well  are  the  functions  of  tem- 
perature and  pressure  (density).  Atypical  diagram  of  the  reduced  viscosity  rj  * in  terms  of  reduced  den- 
sity p * and  temperature  T*  is  shown  on  Figure  5.  Areas  of  moderate  dense  gas  (G)  and  liquid  (L)  are 
pointed  out.  The  line  relating  to  the  diluted  gas  limit  (DGL),  the  saturation  line  (SL)  and  the  melting  line 
(ML)  are  mapped  on  the  diagram  as  the  bounding  lines  of  the  jj  * - p * -T  * surface.  The  picture  clearly 
shows  the  complex  behavior  of  the  viscosity  in  terms  of  density  and  temperature. 

Generally,  analytical  expressions  are  used  for  application  oriented  calculations.  Hence,  accurate  correla- 
tions are  needed.  The  analytical  generalization  of  experimental  data  includes  property  data  in  the  as  wide 
as  possible  range  of  temperature  and  density.  Therefore,  data  near  critical  areas,  saturation  line  and  melt- 
ing are  very  important  in  any  multi-property  analysis,  especially  when  the  correlations  are  to  be  extended 
onto  areas  which  lack  experimental  data.  Additionally,  data  of  viscosity  and  thermal  conductivity  in  the 
dilute  gas  limit  is  needed  as  well. 
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Fig.  5:  Reduced  viscosity  rf  =tj/tj0  in  terms  of  reduced  temperature,  f =T/TC  and  density,  p =p/pc.  G is 
the  gas  state;  L is  the  liquid  state;  CP  is  the  critical  point;  TP  is  the  triple  point;  SL  is  the  saturation  line; 
ML  is  the  melting  line;  DGL  is  the  dilute  gas  limit 

The  investigation  of  viscosity  and  thermal  conductivity  in  the  regions  of  saturation  and  melting  provides  a 
means  for  extension  into  the  range  of  extra-large  pressures.  Results  obtained  from  fundamental  theory 
studies  are  included  as  well  into  the  treatment  of  data.  Viscosity  and  thermal  conductivity  of  the  main 
components  of  interest  (nitrogen,  oxygen,  hydrogen  and  water)  are  reviewed  here.  Experimental  studies  of 
the  transfer  coefficients  for  most  matters  are  primarily  limited  to  temperatures  below  1000  - 1200  K due 
to  high  chemical  activities  of  the  substances.  Hence,  the  most  extensive  measurements  of  viscosity  and 
thermal  conductivity  have  been  performed  for  less  chemical  active  substances  such  as  nitrogen  or  air. 
Generally,  properties  of  the  most  chemical  active  molecule  oxygen  are  less  known.  Table  2 gives  a sum- 
mary of  the  current  boundaries  of  temperature  and  pressure  of  the  experiment  data  for  viscosity  and  ther- 
mal conductivity  for  some  basic  components. 


Viscosity 

Thermal  conductivity 

■ 

Ratified  Gas 
0.1  MPa 

SL 

Dense  Gas 

Gas, 

0.1  MPa 

SL 

Dense  Gas 

ATX 

ATX 

ATX 

Ap.  MPa 

ATX 

ATX 

ATX 

Ap.  MPa 

um 

300-1000 

14-32.976 

to  1000 

to  100 

to  830 

KESH 

to  1000 

to  100 

n2 

400-2000 

65-124 

to  900 

0.1-80 

273-1373 

65-124 

300-430 

to  10 

-126.5 

-300 

to  200 

-300 

0.1-1000 

02 

400-1960 

mm 

300-525 

0.1-80 

300-1000 

55-150 

-300 

to  360 

Air 

400-2000 

70-128 

300-523 

0.1-80 

300-1200 

70-128 

300-1200 

0.1-100 

-300 

to  360 

m* 

-1000 

273.16-643.16 

273-1000 

0.1-80 

-1000  1 

273.16-643.16 

-300 

0.1-100 

Table  2:  Temperature  and  pressure  limits  reached  in  experimental  investigations  of  viscosity  and  thermal 
conductivity 


The  recommended  data  for  viscosity  and  thermal  conductivity,  as  correlated  tables  data  are  included  in  the 
Standard  Reference  Data  (GOST  Standard,  Russia,  NIST  Data  and  others)  [4-6]. 

All  major  theoretical  approaches  in  use  for  computing  viscosity  and  thermal  conductivity  and  diffusion  at 
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high  pressures  are  quasi-empirical.  Due  to  the  specific  behavior  of  the  transport  properties  at  high  densi- 
ties it  is  recommended  to  perform  predictions  based  on  theoretical  approaches  only  when  experimental 
data  is  missing.  Viscosity  and  thermal  conductivity  coefficients  in  the  limiting  case  of  zero  density  may  be 
calculated  applying  Boltzmann  equation  and  Chapman-Enskog  method  [7]  using  potential  functions  to 
account  for  particle  interaction.  Such  data  are  suitable  for  comparison  in  wide  ranges  of  temperature  at 
low  densities.  The  complementary  data  in  the  limit  of  the  saturation  line  and  the  melting  area  may  be  im- 
portant for  correlation  at  high  density. 

Quasi-theoretical  approaches 

One  of  the  wide-spread  methods  for  the  prediction  of  transport  properties  at  high  densities  is  a particular 
method  of  excess  functions.  Usually,  the  excess  of  viscosity  and  thermal  conductivity  is  represented  by  a 
power  low  [8,9].  The  main  suggestion  is  that  the  expression  for  the  viscosity  or  thermal  conductivity  coef- 
ficient consists  of  the  two  terms,  with  one  of  them  a function  of  temperature  and  the  other  a function  of 
density: 

*l(T,p)  = V0(T)+AV  (p)  (1) 

X(T,p)  = X0(T)  + AX(p)  (2) 

The  usual  practice  is  to  express  the  excessive  values  of  A 7]  ( p)  and  AX  ( p)  in  terms  of  a serial  expan- 
sion of  density: 

Artip)  = axp  + a2p2  + a. ^p3  + a^p 4 + a5p5  (3) 

A X(p)  = hyp  + b2p 2 + b^p2  + b4p4  + b5p5  (4) 

Although  some  results  of  the  correlations  (l)-(4)  are  in  good  agreement  with  experiments,  the  method  of 
excess  function  holds  only  within  certain  limits.  In  reality,  the  functions  AT]  (p)  and  AX  ( p)  are  de- 
rived for  either  normal  or  large  temperatures  and  do  not  hold  for  small  temperatures  and  high  densities. 
The  method  of  excess  viscosity  or  thermal  conductivity  yields  reliable  results  only  in  the  region  of  equi- 
distant isotherms.  But  this  equidistance  is  lost  at  high  density  and  at  low  temperatures  near  the  critical 
area.  The  isotherms  of  the  reduced  viscosity  of  hydrogen  are  shown  in  Figure  6.  It  can  be  seen  that  the  er- 
ror near  the  critical  range  may  exceed  40%. 


rc=0.311,  1 03kg/m3 

Fig.  6:  Excess  viscosity  of  Hydrogen  reduced  to  Fig.  7:  Excess  viscosity  of  Nitrogen 

dilute  gas  limit  in  terms  of  reduced  of  density.  near  critical  point 
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The  behavior  of  the  excess  viscosity  of  nitrogen  near  the  critical  area  is  shown  in  Figure  7.  The  enhance- 
ment of  the  viscosity  along  the  critical  isotherm  above  the  regular  part  is  about  37%.  Despite  these  short- 
comings, it  is  obvious  that  the  method  of  excess  functions  is  quite  useful  in  some  cases. 

Generalization  on  the  base  of  Enskog  equation 

There  are  a variety  of  models  proposed  in  the  form  of  reduced  density  and  reduced  temperature.  The  most 
reasonable  equations  for  the  prediction  of  properties  in  wide  areas  of  temperatures  and  pressures  usually 
base  on  modifications  of  Enskog  kinetic  theory  for  hard  spheres.  Among  them  a unique  equation  for  vis- 
cosity has  been  proposed  [10,1 1].  This  single  equation  describes  the  viscosity  of  both  diluted  gases  and 
liquids.  The  fitting  parameters  of  the  equation  have  been  derived  using  the  most  recent  data  of  viscosity  of 
diluted  gases,  of  the  two  phases  on  saturation  line,  near  the  melting  line  and  for  dense  states.  Details  of 
study  are  omitted  here. 

The  quasi-theoretical  equation  of  viscosity  includes  two  terms,  a kinetic  part  and  a collision  part,  as 

^ ~ ^ kin  + ^ col  ’ ^ 

with  pkin  the  kinetic  part  and  pcol  the  collision  part  of  the  viscosity  taking  into  account  the  different  mecha- 
nisms of  impulse  transfer.  These  parts  are  functions  of  reduced  temperature  and  density  (p/  ps)  with  ps 
the  density  of  the  dense  gas  on  the  line  of  melting. 

A series  of  calculations  using  the  proposed  equation  of  viscosity  have  been  performed  and  the  results  are 
compared  with  direct  computer  simulations.  This  comparison  confirmed  the  validity  of  this  quasi- 
theoretical  equation  [10, 1 1].  A set  of  the  results  obtained  is  shown  in  Figures  8-10.  While  Figures  8a  and 
8b  show  the  isobars  and  isochores  of  viscosity  of  H2,  Figures  9 and  10  present  the  isochores  of  the  viscos- 
ity of  N2  and  02,  respectively.  Generally,  the  deviations  of  the  results  are  in  almost  all  cases  less  than  a 
few  percent. 


Temperature,  K 


Fig.  8a:  Viscosity  isobars  of  H2.  Dotted  Fig.  8b:  Viscosity  isochores  of  H2.  Triplet  Point  (T.P.) 

lines  are  placed  near  critical  pressure  and  Critical  Point  (CP)  are  given  on  the  saturation  line 

Model  of  the  effective  diameter  sphere 

Although  the  empirical  correlation  relies  on  direct  experiments  [11]  is  very  useful  and  provides  the  basis 
for  comparison  and  of  reference,  it  cannot  be  applied  for  cases  such  as  the  dense  states  of  H20  or  for  mix- 
tures of  gases  and  liquids.  Furthermore,  a generalization  of  this  correlation  for  other  substances  than  the 
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experimentally  investigated  ones  is  questionable.  Therefore,  the  only  way  of  overcome  these  kind  of  diffi- 
culties is  to  predict  the  transport  properties  on  the  basis  of  kinetic  theory. 


Fig.  9:  Viscosity  isochore  of  N2.  Fig.  1 0:  Viscosity  isochore  of  02. 


A novel  model  for  the  prediction  of  the  viscosity  of  dense  gases  and  liquids  in  a wide  range  of  density  and 
temperature  is  proposed  [12-13],  which  bases  on  specific  characteristics  of  particle  interaction  in  a dense 
environment. 

Both,  the  kinetic  theory  foundation  and  the  generalization  of  the  Enskog  model  for  dense  gases  and  fluids 
make  use  of  the  model  of  a hard  sphere.  They  are  unreliable  for  real  fluids  because  the  actual  interactions 
between  the  particles  are  not  taken  into  account.  However,  if  the  model  of  effective  diameter  of  a sphere  is 
applied  properly  using  effective  potential  of  the  particle  interaction,  the  transport  properties  of  real  fluids 
may  be  predicted  successfully.  The  formalism  of  the  Enskog  hard  sphere  theory  remains  the  same,  but  the 
real  character  of  particle  interaction  is  considered  properly  by  the  effective  diameter  a ~ . 

ceff 

Molecular  dynamics  (MD)  [14-15]  simulations  show  that  the  effective  pair  potential  in  dense  media  is  a 
truncated  function  of  the  distance  between  the  pairs  due  to  a screening  effect  of  neighboring  particles. 
Therefore,  the  MD  calculations  were  carried  out  based  on  the  Lennard-Jones  model  with  either  a cut-off  at 
some  distance:  r ~ (2. 5-3.0)  <y q or  applying  other  modifications  such  as  the  (exp-6-8)  truncated  model. 

Nevertheless,  any  truncation  at  a constant  distance  is  valid  only  for  liquids. 


The  Enskog  equation  of  viscosity  consists  of  two  terms 

/•  m fkm  <aUn  > + fcol  col  > ' 


(6) 


with  the  effective  kinetic  diameter  (the  parameter  of  impulse  transfer  along  the  particle’s  trajectory) 

and  (? the  effective  collision  diameter  (the  parameter  of  impulse  transfer  by  the  particles  collision).  At 

any  temperature  both  parameters  show  systematic  differences.  Similar  results  are  obtained  for  diffusion 
and  thermal  conductivity. 

A new  model  of  the  effective  kinetic  diameter  is  proposed  which  assumes  that  this  effective  kinetic  di- 
ameter is  a function  of  temperature  and  density.  This  model  describes  the  decrease  of  the  free  path  of  par- 
ticles with  increasing  density.  The  interaction  between  particles  is  taken  as  Lennard-Jones  (LJ)  potential 
with  a shielding  function  to  account  for  environmental  effects.  The  cut-off  parameter  is  determined  in 
terms  of  density  and  temperature  and  takes  the  overlapping  of  long-range  attractive  forces  of  closely 
packed  of  molecules  into  account. 
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We  will  describe  the  interaction  of  particles  taking  into  account  the  overlapping  of  distance  fields  due  to 
close  molecular  packing  at  higher  densities  taking  the  effective  potential  in  form  of  a reference  potential 
model  with  an  exponential  screening  function 


r 1 

( 

f \ 

2 ' 

L Vrep  (r ) + £/ aft  (r)  Jexp 

~ 

r 

b 

V 

V cut  J 
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(7) 


with  bcu,  the  cut-off  parameter,  which  accounts  for  screening  interaction  of  the  particles.  The  mean  dis- 
tance between  the  particles  s - n "m  is  taken  as  a screening  parameter.  Thus, 
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with  n the  number  density,  O ^ the  parameter  of  the  potential  function  (LJ  in  this  case).  In  the  case  of  low 
density  the  potential  is  transformed  into  the  reference  potential  function. 


The  results  of  numerical  simulation  using  the  effective  kinetic  diameter  coincide  with  the  results  of  MD 
calculations  for  the  limit  of  high  density  states  as  well  as  with  kinetic  diameter  simulations  of  a dilute  gas 
in  the  limiting  case  of  DGL.  For  dense  gases  the  effective  kinetic  diameter  results  show  excellent  agree- 
ment with  results  determined  from  viscosity  experiments  for  argon,  nitrogen  and  other  gases. 


For  a given  density,  the  effective  kinetic  diameter  is  derived  from  relation 

, 2 . 

,bp  = -xn  , (9) 

with  Q<2'2’,  the  Q -collision  integrals  of  Cowling,  b the  volume  of  hard  spheres  . 


akin  (r  ■*?>* 
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The  model  of  the  effective  kinetic  diameter  uses  no  empirical  parameters  and  can  therefore  be  applied  for 
viscosity  calculations  when  experimental  data  are  missing.  Since  this  approach  accounts  for  density  varia- 
tions, it  is  most  suited  for  property  calculations  with  varying  densities.  The  model  can  be  used  to  predict 
transport  properties  of  dilute  gases,  dense  gases  and  liquids  on  the  basis  of  a single  equation. 

A comparison  of  our  predicted  N2  viscosities  at  100  MPa  with  data  taken  from  the  NIST  is  shown  in  Fig- 
ures 11a  and  b show.  The  agreement  between  the  direct  theory  (present  calculation)  and  the  correlation 
data  (NIST)  is  quite  reasonable.  The  increase  of  the  deviation  with  decreasing  temperature  at  very  low 
temperatures  may  be  related  to  the  LJ  reference  potential  function  whose  parameters  are  not  reasonably 
suited  for  low  temperatures. 


Temperature,  K 


Fig.  1 1-a:  Viscosity  on  N2 : present  calculation  Fig.  1 1-b:  Error  of  viscosity, 

(solid  curve)  and  NIST  data  (circles)  A=[ rjfpresent)- rj(NIST)]/rj(present) 
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Thermal  non-equilibrium  processes 

The  main  processes  in  cryogenic  liquid  rocket  engine  systems  are  droplet  vaporization  and  spray  combus- 
tion including  liquid  jet  atomization,  spray  formation,  vaporization,  multiphase  flow  mixing,  ignition  and 
combustion  in  a high  pressure  high  temperature  environment.  The  time  scales  of  liquid  jet  atomization, 
evaporation,  convection  and  mixing  are  such  that  oxygen  droplets  may  penetrate  into  the  flame  front  and 
enter  the  reactant  which  is  mainly  superheated  steam  for  the  propellant  combination  LOX/GH2. 

Any  evaporation  process  involves  heat  and  mass  transfer  from  the  hot  surrounding  gases  to  the  droplet 
surface  and  vice  versa.  For  an  accurate  prediction  of  these  fluxes  transport  coefficients  and  thermody- 
namic properties  of  both  fluids  are  needed.  However,  investigations  have  shown  that  not  only  conservation 
equations,  equation  of  state  (EQS),  transport  and  thermodynamic  properties  must  be  modeled  accurately, 
but  also  interface  processes  have  to  be  taken  into  account.  Typically,  thermal  equilibrium  conditions  on 
the  droplet  interface  is  assumed  which  may  not  hold  at  very  large  temperature  and  concentration  gradients. 
Therefore,  a model,  which  accounts  for  thermal  non-equilibrium  at  the  droplet  surface,  is  necessary. 

When  the  radius  r of  a spherical  droplets  is  small  and  comparable  with  the  mean  free  path  X of  the  mole- 
cules in  the  surrounding  gas  ( 0.01  < Kn  < 0.2,  Kn=X2r)  the  discrete  molecular  structure  is  considered  in 
a slip-flow  regime.  Foundations  of  the  basic  Boltzmann  equation  are  valid  for  small  Kn  numbers,  and  the 
Navier-Stokes  equations  (the  governing  equations  for  continuum)  may  be  used  for  all  the  entire  regime 
except  the  layer  close  to  the  interface.  This  Knudsen  layer  is  specified  as  collisionless  molecular  flow  re- 
gime. Its  thickness  is  of  the  order  of  the  mean  free  path  of  the  molecules  and  therefore,  the  individual  be- 
havior of  the  molecules  has  to  be  considered.  One  way  of  correct  simulation  of  the  flow  in  the  entire  do- 
main is  to  extrapolate  the  boundary  conditions.  The  boundary  conditions  on  the  solid  surface  under  con- 
tinuum assumptions  are  modified  by  jump  conditions  of  thermodynamic  and  gasdynamic  values.  The  lin- 
ear extrapolation  of  the  real  temperature  field  to  the  surface  gives  the  unknown  boundary  ambient  gas 
temperature  Te  (0)  as  shown  in  Figure  12. 

The  solution  of  Navier-Stokes  energy  equation  with  this  boundary  condition  Te  (0)  coincides  with  the  true 
solution  outside  the  Knudsen  layer.  The  difference  between  the  liquid  surface  temperature  Tin  and  the  ex- 
trapolated gas  temperature  Te  (0)  on  the  surface  ATj  = Tin  - Te  (0)  is  the  temperature  jump.  The  same  as- 
sumption is  taken  for  concentrations  of  evaporated  components.  Accordingly,  the  difference  AY  = Yg  —Yj 
is  the  concentration  jump. 


Figure  12:  Schematic  of  temperature  distribution  around  a droplet 

We  will  consider  here  a model  for  single  drop  evaporation  in  a stagnant  environment  in  an  infinite  vol- 
ume. The  special  problem  of  the  vaporization  processes  of  liquid  oxygen  droplet  (LOX)  in  superheated 
steam  (GH20)  is  investigated.  The  conservation  equations  are  used  for  oxygen  droplet  evaporating  in 
steam  under  the  following  conditions:  droplet  radius  ~ 10  +100fJm,  temperature  ~ 1500+2000  K and  pres- 
sure 1+10  MPa.  Knudsen  number  falls  in  the  range  0.01  < Kn  < 0.15.  Since  there  are  large  gradients  of 
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temperature  and  concentrations  on  the  droplet  surface,  the  boundary  conditions  for  conservation  equations 
must  be  investigated. 

Droplet  evaporation  modeling  of  LOX  in  water  steam  (GH20)  has  some  specific  characters.  The  critical 
parameters  of  H20  are  considerable  higher  as  compared  to  02.  If  surface  temperature  of  the  droplet  is 
lower  than  100  K then  according  to  the  phase-equilibrium  concept,  the  water  vapor  around  the  droplet 
condenses  and  freezes  at  its  surface.  However,  an  ice  layer  on  a L02  droplet  has  never  been  observed  so 
far  in  experiments.  Evidently,  the  processes  in  the  ultimate  boundary  layer  of  L02  and  GH20  may  be  in 
thermal  non-equilibrium  [16-18]. 

Model  of  quasi  steady  - state  vaporization 

If  the  characteristic  time  of  molecular  transfer  is  smaller  than  characteristic  time  of  droplet  evaporation, 
than  the  fields  of  temperature  and  concentration  are  invariant  at  constant  droplet  radius,  and  the  quasi 
steady  - state  vaporization  model  may  be  used.  It  includes  the  following  assumptions: 

• the  temperature  of  droplet  is  constant,  uniform  and  equal  to  the  wet-bulb  temperature,  i.e.  droplet 
heating  is  neglected; 

• the  ambient  gas  has  negligible  solubility  in  the  liquid  and  only  oxygen  diffuses  from  the  surface; 

• the  radial  motion  of  the  liquid  surface  is  assumed  to  be  small; 

• the  pressure  is  equal  to  the  ambient  pressure  everywhere; 

• radiation,  Dufour  and  Soret  effects  are  negligible; 

• density  of  gas  phase  is  constant. 


The  governing  system  of  equations  is  simplified  further  neglecting  diffusion  of  the  gas  into  the  liquid 
droplet,  chemical  reactions  and  viscous  dissipation.  To  determine  the  temperature  and  concentration 
jumps,  Yalamov’s  approach  is  applied  [19-21]. 


The  resulting  simplified  system  of  equations  is  as  follows  [22  - 23] 
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The  boundary  conditions  for  the  gas  phase  are  as  follows: 
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Tf,  7/  are  temperature  and  concentration  of  oxygen  in  the  gas  phase  on  the  droplet  surface;  Tj,  Y\  are 
temperature  and  concentration  of  oxygen  in  the  liquid  on  the  droplet  surface.  pm  is  the  mass  density  of 
the  mixture,  v the  radial  velocity;  Yit  Dn  are  the  mass  fraction  and  the  binary  diffusion  coefficients  for 
each  species,  Sy  is  a term  which  stands  for  different  effects,  for  example,  Dufour,  Soret  effects,  inertia 
source  term  ,etc..;  hm,  h,  are  the  mixture  and  the  specific  enthalpies,  X m is  the  thermal  conductivity  of  the 


mixture;  AH  = h’“t  - h'"'  is  the  specific  enthalpy  of  vaporization  of  species  ; Ky),K[n)  are  coefficients 
of  temperature  and  concentration  jumps,  K(Tn)  ,K'J]  cross  coefficients. 


Furthermore,  we  will  introduce  the  following  the  non  - dimensional  parameters: 
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After  integration  of  the  system  of  governing  equation  and  some  additional  mathematical  simplifications  an 
analytical  solution  of  the  steady  - state  evaporation  problem  is  derived.  The  equations  for  the  temperature 
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and  concentration  distributions 
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can  be  integrated  analytically  for  D = const". 
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For  the  power  dependence  of  D on  temperature,  the  analytical  solution  depends  on  the  power  of  tem- 
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perature.  For  D = D0  — the  equation  for  concentration  distribution  can  be  written  as 
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This  system  of  the  algebraic  equations  for  the  temperature  and  concentration  fields  in  the  neighborhood  of 
the  droplet  is  closed  by  the  expressions  for#,  , Yf,  which  are  determined  from  the  temperature  and  jumps 
conditions. 

Finally  the  system  of  algebraic  equations  for  determination  V , Yf,  Q"  is  the  following 
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with  Qm , @ , the  integral  mass  and  heat  fluxes  and  / the  heat  of  evaporation. 


To  calculate  the  concentration  and  temperature  jumps  according  to  Yalamov  and  co-workers  [21]  we  sim- 
plify the  jumps  coefficients,  as  given  below: 
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with  n0  - molecular  mixture  concentration,  nni  - molecular  concentration  of  fh  component,  n0  = n0I  + n02, 
m,  - mass  of  a molecule;  Du  - coefficient  of  binary  diffusion,  Dn[,]  - coefficient  of  binary  diffusion  first 

approach;  k - Boltzmann  constant;  T0  - temperature  of  a droplet  and  - collision  integrals. 

The  radius  of  evaporating  droplet  is  determined  by  equation 
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Numerical  results 

A series  of  computations  have  been  performed  to  investigate  the  influence  of  various  physical  phenomena 
on  the  distributions  of  the  temperature  and  the  concentration.  We  computed  concentration  and  temperature 
jumps  on  the  surface  of  an  evaporating  single  oxygen  droplet  in  steam  for  initial  droplet  radii  of  10,  25  and 
50  fjm  for  T = 1500+  2500K  and  P = 0.1  and  4 MPa. 

It  is  found  that  the  initial  droplet  size  as  well  as  the  pressure  and  temperature  of  the  surroundings  have  a 
significant  influence  on  the  droplet  behavior  [17].  Due  to  the  lack  of  experimental  data  for  the  binary  sys- 
tem L0X/GH20  we  have  tested  the  model  comparing  our  results  with  experimental  data  [24]  for  the  sys- 
tem LOX/N2,  see  Figure  12.  The  results  are  in  a good  agreement.  The  predicted  lifetimes  of  LOX  droplets 
in  steam  for  various  pressures  and  steam  temperatures  are  shown  in  Figure  13  for  different  initial  droplet 
radii. 


The  temperature  jumps  on  the  boundary  surface  of  droplet  and  surroundings  were  computed  using  differ- 
ent assumptions  about  thermodynamic  and  transport  properties  behavior.  The  physical  properties  in  terms 
of  temperature  and  density  were  taken  in  the  first  case.  In  a second  case,  the  averaged  mean  values  were 
used  [18].  The  corresponding  temperature  jumps  on  the  droplet  surface  are  demonstrated  in  Figures  14a,b 
for  different  droplet  radii  and  two  ambient  pressures  p = 0.1  and  4 MPa  at  various  initial  steam  tempera- 
tures T0.  It  is  evident  from  the  figures  that  the  usage  of  simple  approximations  for  the  thermodynamic  and 
transport  properties  may  cause  significant  errors. 


Fig.  12:  Averaged  Variation  of  (D/Dnf  over  time  for  Fig.  13:  Droplet  lifetime  over  initial  droplet  radius 
the  system  LOX/N2;  P = 0.1  MPa,  7V?  = 300K  for  the  system  LOX/GH2  for  various  pressures  and 

initial  steam  temperatures 

a)  T qh2o =1500  K andb)  Tgh2o=1000  K 
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Fig.  14  a,b:  Temperature  jumps  on  the  boundary  surface  of  droplet  and  surroundings  for  a)  P = 0.1  MPa, 
Tin  — 90  K and  b)  P = 4 MPa,  Tin  = 150  K:  Solid  curves  (1)  are  the  computations  with  real  properties;  dot- 
ted curves  (2)  are  the  computations  with  averaged  constant  properties  at  various  droplet  radii  R.  T0  is  the 
initial  temperature  of  surroundings,  Tin  is  the  temperature  of  liquid  droplet  surface 

Conclusions 

It  has  been  shown  that  viscosity  and  thermal  conductivity  of  fluids  strongly  depend  on  density  and  tem- 
perature and  that  the  prediction  of  transport  properties  which  bases  on  simple  models  of  power  expansion 
of  density  yield  erroneous  results  if  applied  beyond  certain  limits.  The  proposed  effective  kinetic  diameter 
model  makes  it  possible  to  describe  the  viscosity  of  liquids  and  gases  in  wide  ranges  of  temperature  and 
density.  The  performed  numerical  studies  confirmed  that  the  new  model  is  the  less  computing  resources 
compared  with  the  MD  simulations,  which  requires  many  hours  of  computations  per  individual  data  point 
and  many  days  for  mixtures.  Computations,  which  base  on  the  effective  kinetic  diameter  model,  can  be 
done  with  less  effort. 

The  vaporisation  model  which  accounts  for  thermal  non-equilibrium  boundary  conditions  at  the  droplet 
surface  and  describes  thermodynamic  processes  on  the  interface  more  accurately  but  relies  heavily  on  ex- 
act transport  properties. 
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